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' ^ ' Using data from a large-scale three-dimensional simulation of supersonic isothermal tur- 

bulence, we have tested the validity of an exact flux relation derived analytically from 
the Navier-Stokes equation by Falkovich, Fouxon and Oz [2010 New relations for corre- 
lation functions in Navier-Stokes turbulence. J. Fluid Mech. 644^ 465]. That relation, for 
compressible barotropic fluids, was derived assuming turbulence generated by a large- 
c/^ ■ scale force. However, compressible turbulence in simulations is usually initialized and 

maintained by a large-scale acceleration, as in gravity-driven astrophysical flows. We 
present a new approximate flux relation for isothermal turbulence driven by a large-scale 
acceleration, and find it in reasonable agreement with the simulation results. 



1. Introduction 

An important rigorously derived resul t in the statistical theory of incompressible tur- 



bulence is Kolmogorov's four-fifths law (|Kolmogorovlll94lf ) 



{[5u\\{r)r) = -jer, (1.1) 

where 

Su\\{r) = [u{r) - u{0)] • r/r (1.2) 

is the longitudinal velocity difference between two points, and r, separated by a distance 
r = |r|, and e is the mean energy dissipation rate. The four-fifths law is well supported 
experimentally and has been traditionally interpreted as a signa ture of a direct kinetic 
energy cascade within the inertial range of scales. Recent work bv lFalkovich et all ( 201QI 



hereafter FFO), reinterprets the Kolmogorov relation in terms of currents and densities 
of the conserved quantities. From this new perspective, FFO derived an analytic scaling 
relation for barotropic flows with an arbitrary degree of compressibility, based on the 
correlation of currents and conserved fluxes in the inertial range. 

There are few analytic relations proposed for compressible turbul ence, e.g., the rela- 



tion f or an isothermal gas based on an effective energy transfer rate (jGaltier fc Baneriee 



2OIID, and the lack of experimental results at high turbulent Mach numbers makes veri- 



flcation of su ch relations challengin g. While supersonic turbulence is observed in molec- 



ular clouds (jHever fc BruntM2004f ). and is believed to play an important role in star 
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formation (jMcKee fc Ostrikerll2QQ7l ). observations are limited by a vailable angular reso- 
lution, projection and finite optical depth effects (jElmegreen fc Scalo 2004) . At present, 
numer ical experiments provide the best route for testing theories o f compressible turbu- 



lence jS vtine et al' 2000; 'Benzi et al\ l2008l : iKitsionas eraZl 120091 : iKritsuk et ah 12011 



Brandenburg fc Nordlun d 2011 .). Three-dimensional numerical simulations show that 
even at very high Mach numbers, a compressible energy cascade can be recovered with 



200730; IPan et a/.|[2QQ9l ). lAluie et al 



a prop er density weighting, p^l'^u ([Kritsuk et 
( 20121 ) have recently presented evidence from simulations for pressure-dilatation action 
at large scales, and kinetic energy transfer to the dissipation scale via a conservative 
cascade, supporting an earl ier analytica l proof of the locality of kinetic energy transfer 
in compressible turbulence (|Aluiell201l[ ). 

In this article, we use data from a simulation of supersonic isothermal turbulence to an- 
alyze the exact relation for compressible fluids from FFO. We are particularly interested 
in evaluating the role of the properties of the force that acts as a source of momentum 
and energy to compensate for the dissipation losses. We discuss the limitations implied 
by the assumptions used in the FFO derivation and present a new approximate scaling 
relation appropriate for the dri ving commonly used in numerical experiments on com- 
press i ble fluid turbulen ce (e.g., 
20081 : IWang et a/.ll2010h . 



Porter et al 2002; Kritsuk et al 2007a; Schmidt et al 



2. Flux Correlations 

A general relation is derived in FFO, along with a particular case for compressible 
turbulence in a barotropic fluid, both of which are based on two-point statistics in the 
inert ial range. The lower limit of r where the relations are expected to hold is the scale 
where dissipation is negligible, while the upper limit is the correlation length of the force. 
The correlations involve the densities g^, currents j^, fluxes ^ and sources (forcing) /^; 
the sources are assumed to be random, statistically stationary, spatially homogeneous, 
and isotropic. Beginning with the governing equations, and eliminating terms tied to 
dissipation, a relation involving the densities, fluxes, and forcing is derived: 

V, (g^O, mir, t)) = (9^0, t)r (r-, t)) , (2.1) 

see FFO for more detail. 

In order to eliminate the scale dependence on the right-hand side of (j2.ip . the authors 
assume that on scales much smaller than the correlation length (r <C L), the force is 
approximately constant, 

r{0,t) ^ r{r,t). (2.2) 

Therefore, in the inertial range, the correlation of a density and its source is also approx- 
imately constant, 

{q-{0, t)r{r, t)) « (9^0, t)r{0, t)) = la. (2.3) 

Note that this ansatz is essentially carried over from the incompressible case, where 
there is no distinction between the large-scale acceleration and large-scale force since the 
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density is constant Substituting Ca into (j2.ip gives the correlation function 

V, (g<'(0)Ff(r)) = ea, (2.4) 
and assuming isotropy, FFO derive the exact scahng relation in a vector form: 

(q'^mtir)) = (2.5) 

where d is the number of spatial dimensions. 

In the case of a barotropic fluid, the densities are q = (pu^p), and the fluxes are 
F- = puiUj-\-p{p)5ij for z, j = 1, . . . , while Ff~^^ = pUi. When dissipation is neglected, 
the governing equations are the forced Euler equations: 

dtp^V •{pu) = 0, (2.6) 

dtipui) -^djipuiUj -^pSij) = f. (2.7) 

Substituting the densities, fluxes, and forcing into (j2.5p gives a relation for isothermal 
turbulence when p = p, 

Er 

$(r) = ([p(0)w(0) • u{r)p{r)] n||(r)) + {p{0)u\\{0)p{r)) = -, (2.8) 

where u\\ = u ■ r/r is the longitudinal velocity, 

e={p{0)u{0)-f{0)), (2.9) 

and we have taken the scalar product of both sides of (|2.5p with r/r to get (|2.8p in a more 
convenient scalar form. Unlike e in (|l.ip . e is the injection rate of momentum squared, 
not energy. However, in the incompressible limit, the vector form of (j2.8p reduces to a 
third-order velocity correlation function implying (jl.ip . 

In numerical experiments of compressible turbulence, the flow is traditionally driven 
by a large-scale acceleration a, such that / = pa, while the exact flux relation ()2.8p 
was derived under the assumption that / is a smooth, large-scale, field. As noted in 
FFO, physically, / may come from a gradient of potential (e.g., external gravitational 
acceleration, i.e. / = p\/(j) or a = V(/)). The decorrelation and small-scale variations of 
the density at high Mach number make this an important point, as acceleration-driven 
turbulence is important in astrophysics and other areas, such as turbulent convection. 

Let us suggest a new flux relation appropriate for this case. Now, instead of the constant 
e in (|2.9p we shall have the two-point fourth-order correlation function {p{0) p{r)u{0) • 
a(r)), which is scale-dependent. The character of this scale-dependence can be found in 
the particular case of an acceleration short-correlated in time, when {ai{t)aj (0)) = £ij5{t). 
The average of the product of any quantity U {a} and a white noise a is expressed by the 
formula of Gaussian integration via a variational derivative: {U{a}ai) = eij{SU/Saj). 
In the correlation function we consider, it is the velocity field which is related to the 
acceleration (by the equation dui / dt — UiV • u = ai) so that {Sui{t)/Saj{t')) = 5ij0{t — t') 
where the last factor is the step function. We now obtain 

{p{0)p{r)u{Q) ■ a(r)> = {p{^)p{r)) e (2.10) 

where e = (u(0) • a(0)). Since the acceleration is a large-scale field, then a two-point 

f While we recognize that the assumption of lar ge-scale force per unit vol ume is possibly the 
only viable path to a rigorously derived relation (|Galtier Banerje3l201ll . also relied on the 
same constraint to come up with an alternative relation for a compressible case), we emphasize 
that driving variable- density flows with a large-scale force would have a nontrivial effect on the 
turbulence statistics and eliminate the inertial range in the traditional sense, with respect to 
the fluid velocity. 
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velocity-acceleration correlation function can be replaced by a single-point one. Let us 
stress that the constant e is neither the total energy input nor the work of the external 
acceleration. 

Of course, in the physical situations of interest, as well as in simulations, the accelera- 
tion cannot be considered short-correlated in time. In this case, the effective decoupling 
of the fourth moment into the product of second moments expressed by ()2.1Qp is not an 
exact relation but can be suggested as a plausible approximation. We thus come to the 
generalization of the flux relation in the following form (up to an order-unity constant): 

VMr) - (P(0)p(r))£-. (2.11) 

For isothermal turbulence driven by a large-scale acceleration, (|2.1ip provides an ap- 
proximate relation connecting the scaling of ^(r) with that of (p(O)p(r)) in the intertial 
range. 



3. Experimental Verification 

To appraise the relations from FFO, we have used data from a simul ation designed 



to stu dy the inertial range statistics of highly compressible turbulence (jKritsuk et 



20073) • The simulation w as performed with an impleme ntation of the piecew ise parabolic 



method (jColella fc Wood ward 1984) in the Enzo code (|Norman et a/.ll2007h . The forced 



Euler equations for an isothermal gas, (j2.6p and (j2.7p . were solved numerically in a cubic 
periodic domain with a linear size of L = 1 along each axis, and covered by a uniform 
Cartesian grid of 1024^ zones, each having sides of length A = L/1024. The turbulent 
rms Mach number, M = 6, and a steady rate of energy injection were maintained by 
a random acceleration field with power limited to wavenumbers k/kmin ^ [l?^], where 
kmin = 27r/L. The spatially fixed acceleration field was normalized at every time step so 
that {pu • a) was constant; the normalization factor had a standard deviation of ~ 5% 
during the simulation. For this work, we used 43 full data snapshots evenly distributed 
in time in the range t/r G [6,10], where the flow crossing (large eddy turnover) time, 
r = L/2M, assumes the sound speed of unity. For each snapshot, we computed e, £, and 
evaluated <l>(r) for 16 discrete r- values in the range r/A G [8, 128], using 2^^ ^ 4 x 10^ 
random point pairs for each value of r. 

To evaluate the role of various assumptions concerning the forcing and assess the 
hypotheses employed by FFO, we compare $ based on the exact forms (|2.4p and (|2.8p 
with an approximate expression based on (|2.1ip . We begin by examining the spatial 
distribution of the force and acceleration in the lower right-hand panels of figure [H 
which show the values of the x component of a and in a slice through the simulation 
volume at X = 0, where = p/pa is the force normalized by the average density. A 
more quantitative assessment can be found in figure O which shows the power spectra 
of both a and /^ The power spectrum of a falls off at k/kmin > 2, while the spectrum 
of f slowly decreases with /c, and the two are separated by several orders of magnitude 
above k/kmin ^ 2. Further, we measure e = 589, e = 149, (p^) = 3.74, and find that 

= 0.95 e, (3.1) 

supporting the decoupling suggested in (|2.10p . This decoupling, combined with the fea- 
tures of / shown by figures [1] and [21 indicates that the assumption of an approximately 
constant force does not hold in the simulation, as expected. 

Next, we will compare ^ to the scaling from ()2.1ip . We find that the density correlation 
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Figure 1. Slices of several instantaneous dynamical fields through the plane x = 0: (a) logarithm 
of the normalized density field logiQ^p/ p) ; (b) velocity field Ux (velocity is normal to the image 
plane); (c) normalized momentum p'^ = p/pux] (d) acceleration field a^; (e) normalized force 
fx — pI P (^x- N.B.: To account for the distribution of the density field, the colour maps used for 
p'x and fx are highly compressed around 0. 
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Figure 2. Power spectra of the driving acceleration, a (solid black line), and the normalized 
force, f' — pa/p (solid grey line), from a single flow snapshot. 
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Figure 3. Dots: time-averaged density correlation; the standard error of the mean of (yo(O)p(r)) 
is < 0.02. Dashed Une: best flt to the data points with a\og{r/b) in the range r/A 6 [20, 90], as 
indicated by the vertical grey lines. 
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Figure 4. Dots: time- averaged measure of $(r); the error bars show the standard error of the 
mean. Solid line: expected scaling based on a large-scale force. Dashed line: expected approximate 
scaling based on the density correlation. Note that the dashed line is not a fit to the $(r) points, 
but shows the right-hand side of ()3.3|) with parameters determined by the best fit to the density 
correlation ()3.2|) shown in figure [3] 



is well approximated by a logarithmic form: 

(p(0)p(r))^alog(^) (3.2) 

with a = -0.46 ± 0.001 and b = 2.04 ± 0.01 when fitted in the range r/A G [20,90] 
(see figure [3]) . This is expected from the density pow er spectrum approxini ately following 
~ at this particular turbulent Mach number ( Kritsuk et al. l2Q07ah . Substituting 



T2]) into (j2.1ip and integrating radially from to r assuming isotropy, we have 

$(r)c.^[31ogQ-l]r. (3.3) 

Note that ()3.3p is a good approximation only for turbulence at M = 6. Simulations with 
different Mach numbers will have density power spectra with different slopes, and their 
density correlations will have different forms. 

Figure m shows the time- averaged ^(r), the expected scaling for a large-scale force using 
(|2.8p . and the expected scaling based on a large-scale acceleration and our approximation 
(|3.3p . The accuracy of the approximate flux relation can be assessed in figure \5\ where 
we show $(r) normalized by the right-hand side of (j3.3p . using the values fitted to the 
density correlation. The right-hand side of ()3.3p follows the actual measured $(r) within 
±5% in the range r/A G [10,100], somewhat underestimating $ at the lower end of 
the range, and overestimating it as r increases; this behaviour can be predicted as the 
assumed approximations break down for two different reasons. First, from (|3.ip . the 
density correlation and injection rate of specific kinetic energy are not totally decoupled; 
this may account for (|2.1ip underestimating $ at small scales. Second, as r approaches 
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Figure 5. Time- averaged measure of $(r) normalized by the expected approximate 
scaling ()3.3|) . The dashed grey lines show d=5% tolerance levels. 

L/2, (n(0) • a{r)) 0, and as (j2.2p breaks down, a constant e on the right-hand side 
of ()2.1ip causes it to excede $(r). This departure by <l>(r) from ()3.3p at large scales is 
determined by the numerical resolution and driving scales; given a similar simulation 
with higher numerical resolution, we would expect the agreement to scale relative to the 
smallest driving scale. 



4. Conclusion 

We find that the approximate relation we derived for supersonic isothermal turbulence 
driven by a large-scale acceleration field (j2.1ip agrees reasonably well with the simula- 
tion results within a s ufficiently wide range of scales consistent with the inertial range 



presence established in lKritsuk et al\ (|20Q7al ). This implies support for the original FFO 
exact relation (j2.ip . so it can be anticipated that in a system where (j2.3p holds, the exact 
relation ()2.8p will be satisfied as well. Therefore, figures 4 and 5 demonstrate the strong 
sensitivity of the flux correlations in supersonic isothermal turbulence to the nature of 
forcing supporting the statistically stationary state. Such sensitivity has been observed 
in a similar context of the forced one-dimensional Burgers equation, where an invasive 
random force with a oc spectrum ( compare with the grey line i n figure [2]) would pro 



duce a non-universal bifractal scaling (|Mitra et "^ l2005h . Likewise, iBahraminasab et 



(|2QQ8l ) showed that the same Burgers system driven by a force with large-scale correla- 
tion in space and with a Wiener scaling in time exhibits a non-universal scaling for time 
increments larger than the correlation time of the force. Comparison with other relations 
and results addressing energy tran sfer in compressible turbulence ( Galtier fc Banerje3 



20111 : lAluidl201ll : lAluie et a/.ll2012f ) will help us to understand whether non- universality. 



observed here due to the driving mechanism, is a more general property of highly com- 
pressible flows. 
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